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1. INTRODUCTION 

Induction motors (IMs) are used in many industrial applications with wide power range from several 
hundreds to multi-megawatts drives because they are robust and own low cost/power and high power/weight 
ratios [1]. Two strategies of electromagnetic torque control: vector control (VC) and direct torque control 
(DTC) can be used in high precision adjustable-speed IM drives. These strategies give comparable properties. 
Direct torque-controlled drives were introduced more than 10 years later than vector-controlled drives [2]. The 
advantage of the DTC methods is their simple control structures that directly control the torque without many 
frame transformations [3], and therefore, they are simple to simulate on computers and easy to implement on 
practical control systems. DTC strategy brings desired electromagnetic torque control and robustness for 
controlled systems [2]. There are many versions of DTC such as Takahashi’s DTC method (T-DTC) [4-6], 
method with dividing locus of stator flux phasor into twelve sectors (TS-DTC) [7], method that voltage vector 
components are proportional to deviations of flux and torque (DVC-DTC) [8], pulse-width modulation direct 
torque control (PWM-DTC) [2]. 
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The switching frequency is always constant with PWM-DTC. Besides that, the PWM technique also 
ensures the reliable IM excitation and limitation of the stator current vector, reduces switching losses, ripple 
of flux and torque [9]. Current sensors are important parts of DTC drive structures, and processed by suitable 
algorithms. In the paper, components of stator current vector are assumed to be distorded by white noises, and 
in such cases, Kalman filters seem to be the most effective tools. 

Kalman filter (KF) which was invented 60 years ago [10, 11], was applied to smooth noised quantities 
in various fields [12] for example charging state estimation of large-scale battery energy storage systems [13], 
mobile robot navigation [14], impedance parameters estimation for medium transmission line [15], estimation 
of the angle between receiver orientation and receiver-transmitter line in LED communication system [16], 
wildfire progress estimation [17], dimension reduction in X-ray reconstructions of undersampled dynamic 
X-ray tomography system [18]. The Kalman filter brings the optimal estimators for linear systems with 
additive independent Gaussian process and measurement noises. In practice, most control systems are 
nonlinear, extended Kalman filter (EKF) which utilizes Taylor series expansions to linearize a nonlinear 
dynamical model about a working point was made [19]. For highly nonlinear systems, unscented Kalman filter 
(UKF) uses the principle that a set of discretely sampled points can be utilized to parameterize mean and 
covariance without linearization steps [20]. In case of unknown or highly non-Gaussian inputs for linear 
systems, Kitanidis Kalman filter (KKF) was developed [21]. Extended version of Kitanidis Kalman filter 
(EKKF) for nonlinear systems has been also used for state and parameter estimation [22]. State and noise 
covariance matrices in EKF were selected by differential evolution algorithms [23]. In design of stator current 
Kalman filters, discrete-time models of induction generator were obtained by utilizing Euler difference 
method [24, 25]. Flux and speed are estimated using flux models-based EKF [26]. Faults detection are 
implemented by an IM model-based EKF [27]. Genetic algorithm is combined with an IM model-based EKF 
to tune noise matrices [28]. Cubature Kalman Filter estimates load torque using IM state space model [29]. An 
IM model-based adaptive algorithm is added to update system noise covariance matrix in EKF [30]. Speed, 
load torque, and efficiency are estimated by an IM model-based EKF [31]. Stator currents, rotor fluxes, load 
including viscous friction are estimated using EKF and UKF algorithms [32]. In the paper, primitive KF is 
utilized for the filtration of stator current vector components in case of unknown IM model. The paper is 
organized in the section structure: introduction-control structure of induction motor drive using Kalman 
filter-simulation results-conclusions. 


2. CONTROL STRUCTURE OF INDUCTION MOTOR DRIVE USING KALMAN FILTER 
Figure 1 shows the proposed IM drive structure using Kalman Filter. This structure is modified from 

one in [2]: with insertion of KF block between T3/2 block and signal calculation block. The block computes 

orienting angle y, magnitude of stator flux vector w, and electromagnetic torque Te according to (1-5): 


Wsa = J(Usa — Rsisa F )dt (1) 
Wsp = f (usg — Rsisg r )dt (2) 
y = tan" (Wsa/Wsp) (3) 
ps = fe + Wp (4) 
Te = 1.5p(Wsaisg_r — Wsplsar) (5) 


In Figure 1, Kalman filter is used to estimate stator current vector components, and other blocks were described 
in [2]. Because of unknown IM mathematical model, the evolution of these components from time k — 1 to 
time k is simplified according to (6): 


Xk = FXg-1 + Wg-1 (6) 


where x = [ise isgļ": state vector, F is the state transition matrix applied to the previous state vector X,_1, 
Wķ—1 is the noise vector that is assumed to be zero-mean Gaussian with the covariance Q = øĝI. The 
relationship between the state vector and its measurement y at the current time step k is expressed by: 


Ve = HXk + Vk (7) 
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where H is the measurement matrix, Vg 1s the measurement noise vector that is assumed to be zero-mean 
Gaussian with the covariance R = 0,4]. Prediction and update stages of Kalman filter algorithm are computed 


according to (8), (9) and (10-13) respectively: 
Xp = Fagri 
P, = FÊ, FT +Q 
Ze = Yr — HX, 
Ky = PHT (HP, HT +R) 
Xp = Xe + KkŽęk 


P, = (I — Kk H) Pp 


(8) 
(9) 
(10) 
(11) 
(12) 
(13) 


where £: estimated state vector; and Xx: predicted state vector. 
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Figure 1. Control structure of PWM-DTC IM drive with Kalman filter 


3. SIMULATION RESULTS 


Two drive structures without and with Kalman filtering were simulated using Matlab-Simulink. 
Table 1 lists parameters of simulated induction motor. In order to obtain rated voltage, DC link is chosen with 
voltage Vdc = 540 V. Speed controller has proportional gain Kp = 1.5, integral time constant 77 = 0.05 s, limits 
of output +10 Nm. The controller is tuned by finding critical value of the proportional component [33], and the 
integral component of the speed controller is set to zero whenever one of two limits of reference torque is 
reached. The switching frequency of space vector PWM is 20 kHz. Time graphs of physical quantities were 
archieved from the simulated structures at the jump of load torque Jrz (see Figure 2) in case of different values 


of of and of. 


Table 1. Parameters of simulated induction motor 


Parameter 
Rated power 
Rated speed 
Rated voltage 
Rated torque 


Number of pole pairs 


Moment of inertia 
Stator resistance 
Stator inductance 


Magnetizing inductance 


Rotor resistance 


Rotor time constant 


Value 

2.2 kW 

1420 rpm 
230 V/400 V 
14.8 Nm 

2 

0.0047 kg m? 
3.179 

0.209 H 
0.192 H 
2.118 
0.0987s 


The integral time absolute error (ITAE) for two IM drive structures are listed in Tables 2-8. For 
simplicity, abbreviations KF, NF respectively denote drive structures with Kalman filtering and without 
Kalman filtering. In most cases, ITAE for KF is smaller than that for NF. There are 14 cases that ITAE for KF 
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is larger than 75% of that for NF, they all are with small values of of and oå. There are 120 cases that ITAE 
for KF is smaller than a half of that for NF, especially in Table 8, ITAE for KF is only 5.3% of that for NF in 
case of of = oå = 2. Courses of speed and torque, stator current, and stator flux in the case are shown in 


Figures 3-5. 


6p = 0, 
NF 


on =0 0.0349 


o2 =0.25 0.0741 
o2 =0.5 0.1280 
o2 =1.0 0.2386 


o2 =2.0 0.4685 


Gs =0; 
NF 


Oo” =0 0.0475 


o2 =0.25 0.0879 
o2 =05 0.1423 
o2 =1.0 0.2568 


og =2.0 0.4863 


os = 0, 
NF 


KF 
0.0346 
0.0436 
0.0585 
0.0886 
0.1550 


KF 
0.0474 
0.0568 
0.0716 
0.1022 
0.1695 


KF 


Table 2. ITAE in case of Jr, = 0 Nm 


Gs = 0) 67 =0.25 0; =0.25, c0; = 0.5,. 05 = 0:5; -0 = 1.0; or = 1.0, 05 = 2.0; 


NF 
0.0725 
0.1253 
0.1795 
0.2945 
0.5210 


KF 
0.0377 
0.0506 
0.0659 
0.0962 
0.1614 


NF 
0.1239 
0.1778 
0.2344 
0.3488 
0.5723 


KF 
0.0455 
0.0592 
0.0725 
0.1048 
0.1698 


NF 
0.2322 
0.2896 
0.3463 
0.4606 
0.6752 


Table 3. ITAE in case of Jr, = 1 Nm 


65 = 0, Gs = 0.25, 65 = 0:25; of = 0:50, = 05,67 = 1.0, 65 = 1.0; of = 2.0, 


NF 
0.0863 
0.1396 
0.1966 
0.3121 
0.5366 


KF 
0.0513 
0.0639 
0.0793 
0.1099 
0.1756 


NF 
0.1383 
0.1946 
0.2524 
0.3675 
0.5900 


KF 
0.0588 
0.0726 
0.0863 
0.1187 
0.1846 


NF 
0.2501 
0.3069 
0.3644 
0.4781 
0.6984 


Table 4. ITAE in case of Jr, = 3 Nm 


gp = 0. op =0.25, 0p = 0:25; of = 0.5, 65 = 0.5, 07 = 1.0,-05 = 1.0,..02 = 2.0, 


NF 


KF 


NF 


KF 


NF 


KF 
0.0636 
0.0771 
0.0916 
0.1242 
0.1861 


KF 
0.0770 
0.0902 
0.1052 
0.1384 
0.2011 


KF 


NF 
0.4565 
0.5115 
0.5647 
0.6700 
0.8781 


NF 
0.4737 
0.5260 
0.5823 
0.6945 
0.9088 


NF 
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Gs = 2.0, 
KF 
0.1069 
0.1190 
0.1333 
0.1618 
0.2246 


op = 20, 
KF 
0.1201 
0.1323 
0.1472 
0.1761 
0.2389 


gp = 2.0, 
KF 


c =0 0.0762 


oz =0.25 0.1151 
o2 =0.5 0.1688 
o2 =1.0 0.2780 


aż =2.0 0.5042 


öp =l, 
NF 


oh =0 0.1014 


oz =0.25 0.1431 
o2 =0.5 0.2010 
oF =1.0 0.3163 


aż =2.0 0.5496 


or =N; 
NF 


0.0762 
0.0848 
0.1006 
0.1322 
0.1997 


KF 
0.1013 
0.1122 
0.1284 
0.1624 
0.2370 


KF 


0.1127 
0.1654 
0.2187 
0.3324 
0.5571 


0.0786 
0.0918 
0.1080 
0.1391 
0.2051 


0.1635 
0.2163 
0.2728 
0.3861 
0.6089 


0.0860 
0.1002 
0.1139 
0.1465 
0.2127 


0.2698 
0.3260 
0.3820 
0.4934 
0.7166 


Table 5. ITAE in case of Jr, = 5 Nm 


öp = 0 or = 025; of = 0.25; 0; = 05,05 =0.5; 0; = L0 65 = 1.0, 65 = 2.0, 


NF 
0.1410 
0.1978 
0.2551 
0.3719 
0.6058 


KF 
0.1051 
0.1192 
0.1360 
0.1712 
0.2434 


NF 
0.1959 
0.2528 
0.3108 
0.4261 
0.6631 


KF 
0.1130 
0.1278 
0.1434 
0.1803 
0.2526 


NF 
0.3077 
0.3656 
0.4223 
0.5388 
0.7749 


Table 6. ITAE in case of Jr, = 7 Nm 


o = 0; or = 0.25, 165 = 025, gp = 05, 65 =05 a=k o mL om 2.0; 


NF 


KF 


NF 


KF 


NF 


0.1035 
0.1168 
0.1314 
0.1641 
0.2273 


KF 
0.1321 
0.1474 
0.1641 
0.2001 
0.2711 


KF 


0.4881 
0.5437 
0.5980 
0.7090 
0.9465 


NF 
0.5333 
0.5905 
0.6494 
0.7670 
0.9962 


NF 


0.1455 
0.1566 
0.1706 
0.1989 
0.2623 


G5 = 20, 
KF 
0.1787 
0.1922 
0.2087 
0.2420 
0.3139 


gp = 2.0, 
KF 


og =0 0.1280 


o2 =0.25 0.1688 
o2 =0.5 0.2228 
oF =1.0 0.3365 


of =2.0 0.5771 


Gs =, 
NF 


on =0 0.1413 


o2 =0.25 0.1822 
o2 =0.5 0.2360 
oF =1.0 0.3516 


0.1279 
0.1398 
0.1567 
0.1914 
0.2646 


KF 
0.1414 
0.1543 
0.1717 
0.2085 
0.2876 


0.1664 
0.2196 
0.2742 
0.3944 
0.6283 


0.1327 
0.1465 
0.1638 
0.1987 
0.2706 


0.2179 
0.2719 
0.3308 
0.4520 
0.6824 


0.1405 
0.1550 
0.1708 
0.2066 
0.2781 


0.3282 
0.3886 
0.4489 
0.5671 
0.7889 


Table 7. ITAE in case of Jr, = 8 Nm 


6p = 0,. or = 0:25, 05 —0.25,. 65 = 0:5, 6s = 05;.0,-— 10, 05 = 10,65: = 2.0, 


NF 
0.1801 
0.2324 
0.2885 
0.4088 
0.6521 


KF 
0.1464 
0.1607 
0.1788 
0.2161 
0.2942 


NF 
0.2311 
0.2862 
0.3454 
0.4649 
0.7154 


KF 
0.1537 
0.1688 
0.1858 
0.2241 
0.3028 


NF 
0.3424 
0.4016 
0.4604 
0.5772 
0.8690 


0.1601 
0.1740 
0.1898 
0.2249 
0.2928 


KF 
0.1729 
0.1879 
0.2048 
0.2437 
0.3194 


0.5635 
0.6170 
0.6725 
0.7812 
1.0198 


NF 
0.5704 
0.6260 
0.6905 
0.8404 
1.3613 


0.2063 
0.2187 
0.2336 
0.2633 
0.3313 


gr = 2.0, 
KF 
0.2207 
0.2344 
0.2511 
0.2845 
0.3598 


og =2.0 0.5926 
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Table 8. ITAE in case of Jr, = 9 Nm 
g =0; Gs =0; o; = 0.25;.05 = 0.25, of =05;.05 = 05;,-65 — 1.0; of = 1.0,.44 =2.0,. 5 = 2.0, 
NF KF NF KF NF KF NF KF NF KF 
Ou = 0 0.1554 0.1557 0.1944 0.1614 0.2504 0.1692 0.3884 0.1899 0.9126 0.2430 
cå =0.25 0.1983 0.1712 0.2533 0.1783 0.3166 0.1871 0.4811 0.2081 1.3394 0.2600 
cå =0.5 0.2588 0.1918 0.3228 0.1995 0.4024 0.2072 0.5922 0.2296 1.9711 0.2807 
cå =1.0 0.4231 0.2365 0.5152 0.2458 0.6301 0.2558 1.0771 0.2778 3.6842 0.3267 
oå =2.0 1.4037 0.3420 1.9313 0.3513 2.5949 0.3629 4.2353 0.3873 8.4152 0.4448 


N 


In order to evaluate the performance of simulated IM drive structures, the ITAE criterion is utilized. 
The ITAE index has the advantages of producing smaller overshoots and oscillations than the integral of the 
absolute error (IAE) or the integral square error (ISE) indices [34]. In this situation, it is modified according 
to (14): 


ITAE = f’ tle,,(t)|dt (14) 





Reference speed [rpm] 
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Time [s] 
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0.2 0.4 0.6 0.8 l 1.2 1.4 1.6 1.8 
Time [s] 
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Figure 2. Reference speed (upper), and load torque with Jz, = 5 Nm 
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Figure 3. Speeds (upper), and torques with Jr, = 9 Nm, of = 2, 04, = 2 
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Figure 4. Stator currents with NF (upper) and KF with Jr, = 9 Nm, of = 2, oå = 2 
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It is easy to see that starting duration for NF is shorter than that for KF (see Figure 3). The reason for 
this problem is Kalman filter reduces ripple of stator current (see Figure 4) and even ripple of stator flux 
(see Figure 5), and therefore, it lengthens process of reaching rated value of stator flux and limit of stator 
current. In Figure 2, at time of load activation, large negative instantaneous deviation of motor torque and load 
torque makes motor speed decrease quickly for both NF and KF. After this time, for NF, unfiltered stator 
currents continue to give high ripple of motor torque which indirectly lengthens this decrement (see Figure 3). 
Undershoot after load activation of speed responses for NF and KF is 17.2 rpm and 5.1 rpm respectively. 
Ripples of stator flux magnitude in the first and second half of the course (FSHC) for KF are 6.0% and 8.8% 
smaller respectively than those for NF. Motor torque ripple in the FSHC for KF are reduced by 39.9% and 
25.7% compared to those for NF. Especially, ripple of stator current magnitude in the FSHC for KF is 3.3 times 
and 2.6 times lower respectively than those for NF. 


x _ y 
—wy ,NF 
sa H 






Voy ad 


Stator flux [Wb] 


-0.5 





0 02 04 0.6 0.8 l 1.2 1.4 1.6 1.8 2 
Time [s] 


Stator flux [Wb] 


0 02 £04 0.6 08 1 I3 1A 1.6 1.8 2 
Time [s] 


Figure 5. Stator fluxes with NF (upper) and KF with Jz, = 9 Nm, of = 2,074, = 2 


4. CONCLUSIONS 

The PWM-DTC IM drive structure using simplified Kalman filter for stator current filtration in case 
of unknown IM model was presented in the paper. Simulations were carried out with different values of load 
torque jump, noise covariances. The proposed drive structure gave significantly smaller ITAE performance 
index than the conventional drive structure, especially at high levels of noise covariances. The EKF or UKF 
with knowledge of IM mathematical model can be utilized to obtain higher filtration efficiency. Robust control, 
intelligent control or sensorless control techniques can be applied for IM drive with filtered stator current 
components. 
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